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We study the charmonium properties at finite temperature using quenched lattice QCD 



C^- . simulations. Although a simple potential model analysis indicates no bound state at 



T > 1.05T C , our analyses of the spatial correlation between quark and anti-quark and the 



C\) ' spectral function indicate that a bound-state-like structure may survive even above T c 



1. Introduction 

Investigation of the QCD phase transition at finite temperature and density is one 
Qh! of most important subjects in particle and nuclear physics. To create the quark gluon 
plasma phase (QGP), in which quarks and gluons are deconfined and the chiral symmetry 
is restored, the high energy heavy ion collision experiments have been performed and the 



RHIC experiment at BNL is in progress. The J / ip suppression has been considered as one 
of most important signals of formation of QGP phase. However, theoretical understanding 
of this phenomenon is still not sufficient to draw a conclusive scenario. 

We investigate the charmonium properties at finite temperature using lattice QCD 
simulations which enable us to incorporate the nonperturbative effects of QCD from 
the first principle. In this paper, we report our results of three such approaches: (a) 
a nonrelativistic potential model analysis with the static quark potential measured in 
lattice QCD, (6) the spatial correlation between quark and antiquark [I], and (c) direct 
reconstruction of the spectral function from lattice data of correlators using the maximum 
entropy method and the standard x 2 analysis 0. At present stage, all these works are 
done in the quenched approximation (without dynamical quark effects). To circumvent 
severe restriction of degrees of freedom in Euclidean temporal direction, we employ the 
anisotropic lattice on which temporal lattice spacing is finer than the spatial one (for 
setup of anisotropic lattice, see Ref. 0). 

The results of these analyses are summarized as follows. Although the simple potential 
model analysis (a) indicate that there is no bound state above T c , the latter two analyses 
imply that the structure of QCD plasma is much involved and persistent collective modes 
may exist. 
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Figure 1. The static quark potential above T c (left panel), and the critical value fi c and 
measured fi (right). 



2. Potential model analysis 

The static quark potential V^g(r) is determined in lattice QCD simulation. It is ex- 
tracted from the spatial correlation of the Polyakov loops, P(x): 

(P(0)P t (f)) ~ Cexp(-V QQ (?)N t ), P(x) = TrU^U^t). (1) 

As displayed in the left panel of Fig. [3], the results above T c are well fitted to the Yukawa- 
type form, VqQif) = — A exp (— p,\r\) / r , with the screening length 1/fi. The obtained 
value of fi at each temperature is shown in the right panel of Fig. |l|. 

To study whether q-q system has a bound state or not with the obtained potential, we 
solve the Schrodinger equation for a spin averaged state, which is written for the radial 
wave function for the S-wave, $(x), as 

where itlr = m c /2 is the reduced mass and m c = 1.3 GeV is adopted. As the potential 
Vqq(x), we use the Yukawa-type form and obtain the upper bound of the screening mass, 
fi c , below which the bound state is formed (the solid line in Fig. fright)). Since the 
values of fi measured in the lattice simulation are well larger than fi c , we conclude that 
no bound state exists at T > 1.05T C from the the simple potential model analysis with 
the static quark potential obtained in the quenched lattice QCD. This result is consistent 
with Ref. [§. 

3. Spatial q-q correlation 

We examine a spatial q-q correlation (Bethe-Salpeter amplitude) of the charmonium at 
finite temperature |J, 

Mr,t) = EA^ + r,t)Tq{x,t)O\0)), 0(x,t) = (j>(y)q(x + y,t)Fq(x,t), (3) 
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where T = 75 for PS and 7^ for V channels. <fi(y) is a smearing function for which we use 
the result of wr(r, t) at t ^> 1 at T = 0. Figure |] shows the result of q-q correlation in 
the vector channel normalized at the spatial origin, (prij, t) = wr(r,t)/wr(0,t), together 
with the cases of free quark and antiquark. The t-dependence <fir(r,t) give us a hint on 
the behavior of quark and antiquark in plasma phase, since if there is no bound state this 
spatial correlation becomes broader as t increases, like the case of free quarks. Contrary 
to naive expectation, even at T ~ 1.5T C we observe strong spatial correlation similar to 
that below T c , which is stable against t. This implies that the quark and antiquark behave 
differently from the free quarks, and still tend to stay together up to this temperature. 

4. Spectral function 

The spectral function contains the most direct information on the structure of the 
mcsomc correlator, C(t) = Ea?(0(x, t)O t (0, 0)). Recently the maximum entropy method 
(MEM) has been successfully applied to the extraction of the spectral function from 
lattice data at T = without assuming specific form ||. However, its applicability to 
problems at T > is not straightforward, because available physical range of C(t) as 
well as number of data points are severely limited. We propose to use a standard x 2 
fit analysis in combination with MEM, since the former gives more quantitative result 
for the parameters such as mass and width than the latter if presumable form of the 
spectral function is given, which can be estimated with MEM. The numerical simulation 
is performed on lattices with spatial cutoff a" 1 ~ 2 GeV and the anisotropy a a /a T = 4. 

To verify the applicability of these methods at T > 0, we require that they correctly 
work for the correlators at T = with restricted number of data points corresponding 
to the case at T > 0. Unfortunately, MEM does not pass this test for the correlators 
between local operators with present level of statistics (0(500) configurations): MEM fails 
to reproduce the correct feature of the spectral function at low energy even qualitatively 
for the region of C(t) less than O(0. 5fm). The reason may be that one need to extract the 
low energy structure from the correlator in a short distance, where the correlator contains 
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Figure 3. The spectral functions from the x 2 fit analysis at T ~ 1.1T C . 

the contribution from wide range of frequency of the spectral function. This result warns 
against physical significance of MEM analysis of such correlators at T > 0. 

To circumvent this problem, we apply the smearing technique which enhances the low 
energy part of the spectral function. For the correlators smeared with the wave function 
at T = 0, both the procedures works well in the low energy region. Although the features 
of the collective mode such as the mass and width are unchanged by the smearing, the 
smearing can produce a mimic peak structure in the spectral function ||. In order to 
examine this possibility we also observe the correlators smeared with a narrower function 
(half- smeared correlators). The dependence of extracted spectral function on the smearing 
function is compared with that of correlators composed of free quarks. 

At T ~ 0.9T C the results indicate that the structure of the ground state is almost the 
same as at T = 0: almost the same mass with the width consistent with zero. 

At T ~ 1.1T C the result of MEM still exhibits a peak structure around the same energy 
region as at T < T c . Therefore we perform the same kind of the \ 2 fit analysis as at 
T < T c . The fit to the 2-pole ansatz and to the relativistic Breit-Wigner (BW) type 
ansatze give inconsistent results, and the latter fits indicate that the spectral function has 
a peak with almost the same mass as at T < T c and the width of order of 200 MeV, as 
shown in Fig. |3|. The result depends on the smearing function only slightly, in contrast to 
the free quark case, and hence we conclude that this is a physically significant structure. 
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